Conductance distribution in three dimensions: analytic solution of the Generalized 
DMPK equation in the strongly disordered regime 
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We develop a systematic perturbative method to obtain analytic solution of the Generalized 
Dorokhov-Mello-Pereyra-Kumar (DMPK) equation in the strongly disordered regime which de- 
scribes the evolution of the joint probability distribution of the transmission eigenvalues with system 
size. The solution allows us to obtain the distribution of conductance analytically in the insulating 
regime. Our results are consistent with existing numerical simulations of the three dimensional 
tight binding Anderson model, and suggests a possible description of the Anderson transition in the 
presence of a very broad, highly asymmetric conductance distribution. 

PACS numbers: 73.23.-b, 71.30., 72.10. -d 



I. INTRODUCTION 



The apparent existence of metal-insulator transitions 
in two-dimensionsi has generated a lot of interest re- 
cently on the problem of interplay between interaction 
and disorder. However, the consequences of a highly non- 
trivial, non-Gaussian distribution of conductances P{g) 
for even non-interacting electrons at arbitrary disorder 
at zero temperature in two or three dimensions remains 
poorly understood. In particular, several novel features 
in P{g) have recently been discovered2r22, including a 
non-analyticity in P{g) near g = 1 m. the insulating 
regime^, a 'one-sided' log-normal distribution^ near the 
metal-insulator crossover regime in quasi one dimension 
(QID) as well as strong deviations from the expected 
log-normal distribution for an Anderson insulator— li^ 
in three dimensions (3D). They raise serious questions 
about the effect of such broad asymmetric distributions 
on e.g. the Anderson metal-insulator transition in par- 
ticular and on quantum phase transitions"'^- in general. 

Although there is a significant body of numerical work 
availablei^"— , analytic study of the full distribution P{g) 
in 3D at arbitrary disorder is beyond the scope of the con- 
ventional field theory framewor k^T i -.^?, . The novel features 
in P{g) mentioned above were obtained directly from the 
joint probability distribution p(x) of the N transmission 
eigenvalues x = {xq}, Xa > 0, for an A^-channel disor- 
dered conductor, via the Landauer formulaii^ 



P(g)« 



j fldxap{^)6 L -Y,^ec\i\\ . (1.1) 



The N-dimensional integral was evaluated using a saddle 
point method'^''*, while the distribution p(x) was obtained 
from the QID DMPK equation^S and its 3D generaliza- 
tion obtained recently^i2£. For a conductor of fixed cross 
section L^, it describes the evolution with length L^ of 
the joint probability distribution p^^ (x) based on a 'lo- 
cal' maximum entropy ansatz^^. The generalized DMPK 



(GDMPK) equation read: 



,21.22 



N 



9p(x, t) _ 1 v^ _3_ ^ 

1 = 1 



dt 



d dn 

dxi dxi 



pi^,t) (1.2) 



where 



51 = — N 7ij In I sinh Xj — sinh 
i<j 
— \^ In I sinh 2xi \ . 



(1.3) 



In the above, t = L^/l where I is the mean free 
path, jij = 2Kij/Kii where Kij is a phenomenological 
matrix^^ defined in terms of certain eigenvector correla- 
tions that can be explicitly evaluated numerically and 
contains information about dimensionality. The QID 
DMPK equation is recovered when 7^^ — > 1 (we only 
consider orthogonal symmetry). 

Recently a direct numerical solution of Eq. (|1.2I) has 
been obtained^^ by mapping it onto a Langevin equation. 
The numerical results confirm that the GDMPK equation 
indeed distinguishes between two and three dimensional 
models with the same number of transmission channels. 
With detailed comparisons, the study confirms that: (i) 
Eq. (|1.2p correctly describes the full P{g) not only in the 
strongly localized regime but also near the critical regime 
in 3D, and (ii) two parameters, Kn and ^^"12, are enough 
to model the entire matrix Kab at least in these regimes. 
In Rcf. [21] it was shown that at least two parameters are 
needed in order to guarantee the conservation of proba- 
bility, so this can be considered as a minimalistic model. 
Within this approximation, disorder is characterized by 
r = l/KiiLz. In the insulating regime one can interpret 
^ = 4^/1^11 as the localization length, giving P = ^/ALz. 
The other independent parameter 712 = C/8i, so that 
r/712 = "^L/Lz depends only on the geometry of the 
system. Therefore for a cubic system, within this ap- 
proximation, the entire distribution is characterized by a 
single disorder parameter P. 

As shown in Refs. [ll.J2i the parameter Kii{L), where 
L is the cross-sectional dimension, contains informa- 
tion about the dimensionality of the system via the 



dimension-dependent eigenvector correlations that define 
the matrix K . In particular in 3D, 



Kii ex 






insulating regime, 
at critical point, 
metallic regime. 



(1.4) 



such that the quantity Kn = liniL^oo Kii{L) is zero in 
the metallic regime as well as at the critical point, but 
it is finite for insulators. Clearly the parameter Kn can 
be considered as an order parameter for the Anderson 
transition. We can therefore expect that Eq. (11.21) can 
be used as the starting point for studying the Ander- 
son transition in terms of the full distribution of con- 
ductances, starting from the insulating limit. However, 
unlike the QID DMPK equation, a general analytic so- 
lution of the GDMPK equation for arbitrary disorder is 
not yet available. 

In this paper we will obtain an analytic perturbative 
solution of the GDMPK equation for N transmission 
channels valid in the strongly disordered regime F ^ 1. 
The solution is based on a mapping of the equation on 
to an imaginary time Schrodinger equation^^, describing 
N interacting bosons in one dimension evolving from a 
delta function initial condition and satisfying some sin- 
gular boundary conditions. In Ref. |T2|, the full P{g) in 
3D was evaluated based on a solution of the GDMPK 
equation in the limit where the interaction between the 
eigenvalues was totally neglected. In this work, we in- 
clude the interaction and show that it changes the quali- 
tative features of the distribution. We begin by develop- 
ing a systematic perturbation theory method to obtain 
the A^-particle Greens function. As a first check, we show 
that our solutions agree where exact solutions are avail- 
able, namely for the special cases of 712 = 1,2,42^. We 
then use the A^-particle Greens function to obtain the 
joint probability distribution of the eigenvalues Xa > 0, 
a = 1,2 ■ ■ ■ N, analytically. We show that in particular, 
the density of the eigenvalues remains uniform, indepen- 
dent of disorder, with an exponential gap at the origin 
that increases with increasing disorder. Indeed, it is this 
feature which is at the heart of many novel properties in 
the distribution of conductances in 3D. It suggests that 
the Anderson metal to insulator transition can be viewed 
as the opening of a finite gap in the density of the eigen- 
values at the origin. Finally, we do the A'^-dimensional 
integral of Eq. (jl.ip approximately to show that the dis- 
tribution of conductance agrees with the numerical re- 
sults, improving on the approximate solution obtained in 
Ref. 12] where the distribution was found to be much 
broader. In particular, we obtain the variance a{lng) as 
well as the skewness x(ln5). We find that asymptotically 
the variance is not a simple power law although numerical 
data can be fitted with a power law in limited regimes. 
We also find that asymptotically the skewness goes as 
a positive constant. Further, the skewness changes sign 
from positive to negative well before the metal-insulator 
transition point. A short version containing some of the 



above results was published in Ref. [26} . Here we provide 
details of the calculation and also these new results. 

The paper is organized as follows. After a brief in- 
troduction in section I, we discuss the mapping of the 
GDMPK equation on to a Schrodinger Equation in sec- 
tion II. In section HI we show how a systematic perturba- 
tion theory can be developed and point out how a resum- 
mation of the perturbative expansion is essential in order 
to satisfy certain boundary conditions. We will use the 
results obtained in section HI to obtain the joint proba- 
bility distribution of transmission eigenvalues in section 
IV and then the full distribution of conductances in sec- 
tion V. Section VI contains summary and conclusions. In 
Appendix A we discuss how to include systematic correc- 
tions to the leading order results of section HI. 



II. MAPPING GDMPK ON TO A 
SCHRODINGER EQUATION 

Within the approximation that only two independent 
matrix elements Kn and K12 can describe the entire K- 
niatrix, Eq. (II. 2p can be written in the form 



dp{x, t) 
dt 


K -^ d 
4 fr^Sx, 


' d Off 

dxi dxi 


P{ 


n = 


— 7 S ' In sinh^ Xj — sinh 

— SJln sinh2a;i . 



p(x,i) (2.1) 



(2.2) 



where K = Kn, and 7 — 712. The initial condition is 
given by: 



p(x,< = 0) = S{x 
and the boundary conditions are 

" d dn 
dxi dxi 



lim 



p = 0; lim p = 0. 



(2.3) 



(2.4) 



Following Ref. [25[ we map the equation for p on to an 
imaginary time Schrodinger equation. First we use a fac- 
torization 



p(x,t|y)=e(x)GAr(x;t|y)r'(y); C 



-0/2 



(2.5) 



The TV-particle Greens function GAr(x;i|y) = 
G{xi,X2,- ■ ■ ,xj^]t\yi,y2,- ■ ■ tUn) then satisfies the 
imaginary time Schrodinger equation in the variables x^, 
given by 



dGN/dt = HG 



N 



(2.6) 



where 



' 4 dx. 



2 -f u{xi) 



Y,^{x,,x,j). (2.7) 



i<j 



and the single particle potential and the interaction terms 
are given by 



III. THE GREENS FUNCTION 



v{Xi,Xj) 



A,, 



sinh 2.T 



sinh (xi 
1 



Xj) 



sinh {xi 



Here 



A,; 



-K/2 

K"f['-f - 



2)/8 



(2.8) 



(2.9) 



The Greens function G^i'^] i|y) will then have the ini- 
tial condition 

1 ^ 

G^(x; t = 0|y) = ^ J] e5,As(y) J] ^(^^ " 2/0 (2-10) 



^(y) 



where the subscripts S and AS refer to the symmet- 
ric and antisymmetric initial conditions, respectively, de- 
pending on whether fl is symmetric or antisymmetric in 
its variables and n stands for permutations. In the small 
y limit, this recovers the proper initial conditions on p. 
The boundary conditions become 

lim [^ _ -%-] ^ 0. (2.11) 

xi^o [ oxi sinh2a;ij 

Note that the factorization requires defining a new set of 
variables yi, but in the end we will be interested in the 
y — >■ limit 

p(x,i)=p(x,<|y = 0). (2.12) 

Therefore, we will be interested in the y — > limit of the 
iV-particle Greens function Gn{^, t\y) as well. According 
to the construction of Eq. (|2.5p . p(x,t|y) looks highly 
singular in the limit y — > 0, given by 



p(x,i|y) 



Yli^j I sinh^ Xj - sinh^ 2;^! 3 Hi I sinh2xi 



I\i<j I sinh^ Vj - sinh^ yi| 2 J]^ | sinh 2y^\i 

GAr(x,i|y). 



(2.13) 

However, it turns out that the complete A^-particle 
Greens function can be obtained exactly for the special 
values 7 = 1,2,4. In these special cases^^ the small ym 
limit of Gat has the form 



Gjv(x,i|y) 



JJlsinh^j/j - sinh^yi|2 
[||sinh22/,|5; y„<l. (2.14) 



Thus the Gjv for these 7 cancel the singularity arising 
from £,~^{y) exactly as j/m — )■ and yields a finite result 
for p(x, t). This is a highly non-trivial result, and typi- 
cally an approximate evaluation of the A^- particle Greens 
function will not have the same property, rendering the 
extraction of p(x, t) impossible. In what follows, we will 
organize the perturbative series for Gjv in such a fashion 
as to conform to this symmetry. 



As given in Eq. (|2.6p . the A^-particle Greens function 
satisfies a differential equation which is different from 
the usual many-body Greens function defined in terms 
of time-ordered annihilation and creation operators. In 
fact it satisfies the differential equation obeyed by N- 
particle propagator. Taking a cue from this, one can 
straightforwardly verify that Gn can be written as: 

N N 

G^(x,t|y) = — (0| n V'(x™)e-«* [] V'^(yn)|0)(3.1) 



where |0) is the A^-particle vacuum and 
j_Jx^\x)^-- — 

1 nOO nOO 

- dx dx'ip\x)ij\x') 



H 



u{x) 



■0(a;) 



(3.2) 



X v(x,x')ip{x')^j{x). 



This marks a distinction between the Gat of Eq. (|3.ip 
and more typical many-body Green's functions in that 
the expectation in this case is against the vacuum, rather 
than the many-body ground state. As we will see, this 
will have dramatic consequences for the diagrammatic 
expansion of Gjv- Now as usual we would like to go to 
the interaction picture, and so we define the interaction 
picture imaginary time evolution operator: 



S{t2,ti) = Texp 



dt' H'{t') 



(3.3) 



where H' is the perturbation in the interaction picture. 
Taking advantage of the fact that H' annihilates the vac- 
uum, and that t is always greater than 0, we can write 
Eq. dSU as: 

1 ^ 

Gjv(x,i|y) = — (0|T5(oo,0)[] V(a:™,i) 



N 



771=1 



(3.4) 



From the standard theory of non-equilibrium Green's 
functions, the diagrammatic expansion of Gat follows. 
We would compute all diagrams with N external legs Xi, 
and N external legs yi as illustrated in Fig. [TJ And we 
would connect them via the external potential and inter- 
action vertices shown in Fig. [51 Finally, the unperturbed 
Green's function would be given by: 



Go(ii,t2) = -(o|r^(ii)^Hi2)|o) 



(3.5) 



We would integrate these diagrams over all space and 
over time from to 00. But we will observe that Gq is 
zero for negative times ii — ^2, and so this will allow us 
to integrate over all time as well without error. 




XN-l,t 



yN-1,0 



XN,t / \yN,0 

FIG. 1: Diagrammatic expansion of Gat 




-h(i) 



-v(xXlo(t-fl 



FIG. 4: Typical diagram, W, in Taylor expansion of Gn 



greater than 0. 



FIG. 2: Vertices appearing in diagrammatic expansion of Gn 



A. Simplification of the diagrammatic expansion 

The fact that the expectation is taken against the vac- 
uum has a few important consequences for the diagram- 
matic expansion. First we will find that we have no 
self-energy corrections to Go- This is because the self- 
energy is due to the interaction of the particle with the 
background state. But since there are no particles in our 
background state, the self-energy vanishes. Also, because 
Go does not allow negative times, all diagrams must pro- 
ceed forward in time. This will rule out all crossed dia- 
grams, vertex corrections, etc. The vacuum expectation 
dramatically simplifies the diagrammatic expansion, and 
we can determine its general form from the following ar- 
gument. Taking A'^ = 3 for illustrative purposes (the 
argument for general N is completely analogous) , G3 can 
be constructed via the expansion in Fig. [3] 

The proof is straightforward. Simply write out the 
expansion algebraically, and then act on both sides of 
the equation with the operator —d/dt— Hq, where Hq is 
the kinetic energy term. The result will be the equation 



d_ 
dt 



iJo - ^ u{Xi) - ^ V{X^, Xj) 



i<j 



Gs 



S{t)S{xi,X2,X3\yi,y2,y3) (3.6) 
which is indeed the desired differential equation for times 



■ + + — < — + ^' — + — ^ + I 1- 



B. Exponential expansion 

We will find it most convenient to analyze the series 
expansion of the Green's function in terms of an expo- 
nential series (defined below), rather than a Taylor series. 
One reason is that we will find the individual diagrams 
in the Taylor expansion diverge logarithmically (for our 
potentials) for small y. Yet when exponentiated, the di- 
agrams will give logical, finite results. Additionally, we 
will need to extract from the Taylor expansion the non- 
perturbative symmetry required in Eq. (I2.14p . and this is 
perhaps most directly accomplished using the exponen- 
tial series. To that end it is useful to consider the follow- 
ing. Each diagram in the Taylor expansion of Gn can be 
labelled according to the number of single particle poten- 
tial lines on each particle, and the number of interaction 
lines between each pair of particles. For convenience we 
will then label each diagram, W, in the Taylor expansion 
according to the following notation; 



W, 






The superscripts, {ui}, label the number of external po- 
tential lines on each particle, i = 1...N. The subscripts 
{vj} label the number of interaction lines between each 
pair of particles j = l...N{N — l)/2. For instance, the 
diagram in Fig 0] would be labelled WIq2 since it has two 
u's on particle 1, one u on particle 2, and one u on parti- 
cle 3, as well as one v between the first particle pair (12), 
v's between the second particle pair (13), and two v's 
between the third particle pair (23). Other diagrams ob- 
tained by rearranging the potential lines on this diagram 
would be of the same order. Using this notation, the 
Taylor series expansion of the N-particle Green's func- 
tion could be written as 



FIG. 3: Dyson expansion of Gg. Gg is represented by the 
bold lines. 
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(3.7) 



where G)sf is the unperturbed N-particle Green's func- 
tion. In preparation for the next step, we'll divide 

through by G)^' to get 



G 



N 



G 



= 1 



E E 

{Ui}=0{vj} = 






(3.8) 



(0) 



where we define the dimensionless diagram w — W/Gjy 
Now we'd like to re-express this Taylor series as an ex- 
ponential series. So we write, 




FIG. 5: Typical diagram, W, in Taylor expansion of Gs 



G 



N 



G 



= cxp 



E E hvh 



{tii}=0{iij}0 



-< — r 



(3.9) 



To define / with respect to w we equate both expressions 
and solve for / in terms of w order by order, where 'order' 
is understood in the sense of the subscripts/superscripts 
"1 '■■■'"" . This process in general defines f in the 

exponent. Formal rules can be worked out and they are 
as follows: 

The N-particle Green's function is given by: 



Gat = G^^ exp 



E E /. 

{ui}=0{t)j}=0 



{n^} 



(3.10) 



where the superscripts {ui] label the number of exter- 
nal potential lines on each particle with i = 1...N, and 
the subscripts {vj} label the number of interaction lines 
between each pair of particles with j = l...N{N — l)/2. 
There is a one-to-one correspondance between a particu- 
lar dimensionless diagram, w, in the Taylor series expan- 
sion of Gat, and a particular term / in the exponential 
series expansion. The rules for constructing the / corre- 
sponding to the w are as follows: 

Start with the basic diagram w,^ 

Vertically cut w) '( inbetween interaction and po- 
tential lines, into all distinct diagram factorizations 
(daughter diagrams) such that the sum of the cor- 
responding indices in each daughter diagram adds 



1, 



up to: 



3. For each daughter diagram, each w receives a (-) 
sign. 

4. For each daughter diagram, let n be the number of 
if's in the daughter diagram. Then divide by n. 

To illustrate, consider the particular diagram in the 
Taylor series expansion in Fig. ([5]) (for N = 3). We can 
cut this diagram in 3 different ways. One way is between 
the two interaction lines in Fig. ([5]), between the last 
interaction line and the external potential line in Fig. ([T]), 
and between both, Fig. ([5]), 



"OlO "oio - ^ ^^ 



2' X 2^*- 



-^ — 2' 



FIG. 6: Daughter diagram 1 



resulting in the following expression for /: 



1 



I 



/020 ~ ^020 + o ^010^010 + 0^020^000 



- ^W, 



010 000 000 

ooo'"oio"'oio- 



(3.11) 



We recall that w = W/G3 in this case. 



Extraction of symmetry-satsifying component 
of Giv 



There is a set of /'s appearing in the exponential series 
of Gat that we can sum exactly. These are the ones which 
give us the exact single particle Green's function, Gi , and 
two-particle Green's function G2. It turns out that these 
terms are also responsable for the asymptotic behavior 
in Eq. (12.141) . Going back to Eq. p. 101) . consider the set 
of diagrams parameterized by /(^"^ o' °- This set has ui 
external potential lines on particle one, and nothing else 
on the other N-1 particle lines. If we add these diagrams 
up, then we get: 



^^P^ 2^ 2^ /00...0 



= 1 



,,DiO...O 



E E '"00.. .0 

{v.}=0{uj}=0 

Gi(l)G^^i 






Gi(l) 
G?(l)- 



(3.12) 



Similarly, adding up /oa^.'o'° would give us Gi(2)/G5. 
Now consider the set of diagrams given by fy%"'^Q- This 
would correspond to the set of interaction ladders be- 
tween the first pair of particles, say particles 1 and 2. 




FIG. 7: Daughter diagram 2 
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FIG. 8: Daughter diagram 3 



And there would be no interaction lines between any 
other pair of particles (and no external potential lines 
for that matter). If we add these up, we will obtain 



exp • 



E E 

{v,}=0{«,}=0 



(•00. ..0 
JviO...O 



1+ E 

{«j}=0 {ui}=0 

6*2(1, 2)G5/_2 



00 
V^ ,00. ..0 



Z^ ■Wv,0...0 






^2(1, 2) 

G§(1,2)' 



(3.13) 



be defined via Hq + 5u, where 5u is the single particle 
perturbation. G2 would be defined via Hq + v, where v 
is the two-particle interaction. For instance, if the single 
particle potential is simple enough, then the exact single 
particle basis may be convenient to use. In that case, 
there would be no external potential lines appearing in 
the the /' diagrams, only interaction lines. 

It will be evident that the terms to the left of the ex- 
ponential in Eq. p.l4p possess the symmetry present in 
Eq. (|2.14p . and as such would constitute a natural 'first 
order' approximation, though they contain contributions 
up to infinite order in both A„ and A^,. The /' terms in 
the exponent are corrections to this approximation. As 
evidence of the adequacy of this first order approxima- 
tion, we adduce the fact that by itself, it reproduces the 
solution of the QID DMPK equation for the distinctly 
non-perturbative values of 7 = 1,2,4. This in itself is 
remarkable, since the DMPK equation is a many-body 
problem, but the solution, at least in the metallic and in- 
sulating regimes, can apparently be factored into 1 and 
2 body Green's functions. Thus, the first order approx- 
imation may be more generally accurate than it might 
seem. Indeed, the first order approximation treats all 
one and two particle correlations exactly to infinite or- 
der, and approximates all higher order correlations as 
products of these lower order correlations. The /' dia- 
grams correct this approximation. So if three particle 
correlations are important then the /' terms will need 
to be included. But if not, then we may expect that our 
'first' order approximation may be highly accurate - even 
for non-perturbative interactions. 



and similarly for foy^'. o^ ^^c. Therefore let us separate 
the /-terms with only one non-zero index from the rest, 
which we'll denote by /'. After adding up the one-index 
/'s we will obtain. 



Gn = GJ^^exp 



Giv°^n 



E E hvh 

{vj}=0 {ui}=0 

GiWyr^aM 
G?(*)f4G0(z,j) 



X exp 



{vj}=0{uj}=0 



(3.14) 



This is our main result for the Green's function, G]\[. 
Again, /' consists of all diagrams which don't consist 
purely of external potential lines on one particle, or in- 
teractions lines between one pair of particles. The first 
order contribution to /' will be considered in the Ap- 
pendix. 

On a general note, the unperturbed basis, Hq, does not 
have to be the free particle basis; any basis will suffice. 
Given a choice for Ho, Hq will define G}^ , and will be the 
basis in which the diagrams in /' are expanded. Gi would 



IV. DISTRIBUTION OF THE EIGENVALUES 

Now we would like to apply this formalism to our prob- 
lem of calculating p(x, i|y). To wit, we need to calculate 
the 1 and 2 particle Green's functions. Since the sin- 
gle particle potential u{x) is not small, we would like to 
include it in our basis Hamiltonian, Hq. 



A. Single particle Green's function 



Gi will satisfy the differential equation. 



9 M°) 

Ft' 



K d^ 



g(°) 



(4.1) 



This equation has been solved in the literature^. 
The un- normalized eigenfunctions are ipkix) ~ 
sinh xP" (cosh x) where ly = —1/2 -I- ik, and P^{x) 
is a Legendre function. The small/large x behavior of 
the eigenfunctions is given by 



r.1/2 



X << 1 



■>Pk{x) ^ 

ipk{x) ^ VxanhxRe 



r{ik)e""' 



r(ifc + 7/2) 



X > 1. (4.2) 



For large times (lengths) t = Lz/£, only the lowest ly- 
ing eigenstates {k « 1) are important to the Green's 
function, and we will make this approximation on the F 
function so that F(jfc)/F(l/2 + ik) « 1/ik. Normalizing 
the eigenfunction and evaluating the Green's function we 
obtain, in the large x/small y limit: 



N 



G%{^,t\y«l)=l[ 



L^ VKirt 



er^l^'xd,l\ 



(4.3) 



effective boundary conditions imposed by u{x). To sat- 
isfy these, we can add to this solution Green's functions 
corresponding to 'image charges' outside the boundary 
2^1, 2:2 = 0. Then we will have: 

G2(x,f|y) 
= 5(2:1 - 2:2, t|yi - yi)g{xx + xi,t\yx + yi) 

- 5(2:1 -f 2:2, t|yi - yi)g{xx - 2:2, %i + 2/2) 

- g{-xx - X2, t\yi - y2)g{-xi + 2:2, t|yi -I- ^2) 

+ g{-xi + 2:2, t\yi - y2)g{xi - X2,t\yi + ^2). (4.9) 



B. Two-particle Green's function 



Finally, forming the ratio G2/G2 and taking the small 
y limit, we obtain, 



Next we need the two-particle Green's function. Since 02 (x, t\y << 1) 

we're using the exact single particle basis, the external G°(x t\y « 1) 

potential u{x) will appear in the equation for G2, 



Ty{xi + 2;2)r^(|xi - 2:2!) 



X lyl-yl^/' 



(4.10) 



d_ 
dt 



G2 



E 






u{x.i) ) +v{xi,X2) G2 

(4.4) 

This equation would be rather difficult to solve, so we 
will use an approximation afforded to us in our region of 
interest. Note above in Eq. (|4.2l) . that ipk{x) ~ sin(fc2:) 
for 2: > 1 and k << 1. Since we're concerned with the 
insulating state behavior, we will naturally be consid- 
ering only large x. Therefore we may replace the ex- 
ternal potential u{x) with the boundary condition that 
''pkix) = sin(A:2:). Given this, we can evaluate the Hq two- 
particle Green's function in the large x, small y limit. 



Gf(x,i|y«l) 



2xiX2yiy2 
^exp 



Tr{Kt) 
G2 will then be defined through 



Kt 



(4.5) 



where T^{x) is defined as: 



T^(x) 



V2 
and we have defined 



1 ^— erfcp 



2V2 V2Kt 



erfcp(2;) = -^^e"^ erfc(2;) 



(4.11) 



(4.12) 



Left out of Eq. (|4.10p are the additional terms coming 
from the effect of the boundaries. These also have the 
appropriate symmetry, but are negligible when 2; > 1. 
Ultimately, we could have ignored the external potential, 
as we did in our previous paper j^, for purposes of calcu- 
lating G2, since it has effect only near the origin, but in 
the insulating regime, only 2; >> is of consequence. 






K d^ K d^ 



G2 (4.6) 



C. N-particle Green's function and p(x, i) 



such that we use sin(fc2;) as the free basis. This equation 
can now be solved. We first change variables to zi = 
xi — X2 and Z2 =2:1+2:2. Then the equation is separable, 
and may write G2 — g{zi)g{z2) where g{z) satisfes. 



d 

81' = 



K 92 if 7(7-2) 1 

"T9^^T 2 



sinh z 



9- (4.7) 



Analogous to the single particle case, un-normalized 
eigenfunctions are given by 



tpkiz) 



,7/2 



^pkiz) ~ vtanhzRe 



F(ifc) 



z << 1 



F(ifc + 7/2) 



z > 1(4.8) 



To evaluate g{z) in the insulating limit, we let k << 
1 and 7 << 1 as well. With this approximation, 
T(ik)/T{ik + "f/2) w (ik + "f/2)/ik, and upon normaliza- 
tion of the ■0, we can construct the green's function g{z). 
Now the Green's function g{zi)g(z2) does not satisfy the 



From Gi and G2 and using Eq. (|3.14p we can construct 
the approximate N-particle Green's function to first or- 
der. 



^ 1 
G.^U-T^ 



=e-</^*x,x 



Y[T^ix, + x,)T^i\x, - x,\)yl%j - y]yl\ (4.13) 



i<] 



Observe that we have the requisite symmetry in the small 
y limit, as required by Eq. (J2.14I) . Now from this we may 
construct p(x, t) using Eq. (|2.13p . We obtain. 



^ 1 



-.e -"I'^'x. 



-^ i V^sinh 22:i X 

, \JK-nt 

I I I sinh Xj — sinh Xi\^' T^{xi + Xj)T~^{\xi — Xj\). 



(4.14) 



D. Deeply insulating limit of the smallest 
eigenvalue 

We can write the eigenvalue distribution in the form 
of a classical partition function. First recall that the 
'time' variable is the longitudinal length of the conductor 
t = Lz/L So we have, 



1 

'kI 



AIL., 



Now we may write p(x) as, 



p(x) ex exp 



-Y,u{x,)^Y.'^{, 



i<j 



(4.15) 



(4.16) 



where we have defined: 



U{xi) = Fxj —h\Xi lnsinh2xi 



V \Xi , Xj 



-— ■ In I sinh Xi 



sinh Xo 



- \n[T^{xi + Xj)T^{\xi 



(4.17) 



In the insulating limit, 7 << 1, the XiS are large, and 
the conductance is dominated by the smallest eigen- 
value. Therefore the conductance is given hy g ~ 



X;i/cosh" 



1/cosh xi. In this case, we can obtain 



the full distribution P{g) by considering the contribution 
from xi alone. We can approximately obtain P{xi) via 
the following procedure. We order the eigenvalues from 
least to greatest, separate out xi, and then integrate over 
the rest. 



p{xi) ex e-^("i)ZAr_i 



(4.18) 



where 



00 00 00 



Zn-i ^ dx2 dx3... dxNPN-i{x) (4.19) 



xi Xi X\ 



and 



PAr-i(x) = exp[- y^ [U{xi) + V{xi, 

1=2 

- ^ V{x^,Xj)\. 

2<i<] 



(4.20) 



We can consider Z^si-i to be the partition function of a 
classical gas of A^ — 1 particles constrained between xi 
and cxD. Calculating this classical partition function is 
non-trivial in itself. Though the interaction is weak, a 
cluster expansion will not work because the interaction 
is long-ranged; the individual terms in the expansion will 
diverge with the particle number. Barring a more accu- 
rate method, we will use the mean-field approximation^!. 
The basic idea is that 



^Af-l 



F„ 



(^)p-(lnp-^), 



(4.21) 



where {E)p is the expectation of the energy of these par- 
ticles with respect to p, (Inp^^)p is the entropy, and p is 
the probability distribution function for the A^ — 1 par- 
ticles that minimizes the free energy Fp. The mean- field 
approximation consists of postulating a trial p: 



N 



p{x2,X3,...,xn) = Yi 



n{xi) 
N - 1 



(4.22) 



where the n{x) is normalized to N — 1, and then mini- 
mizing Fp with respect to n. Forming Fp with our trial 
p, we obtain, 

00 

Fp = dxn{x)[U{x)+V{xi,x)] 

xi 

00 00 



+ 2 / ^^ / d'yn{x)n{y)V{x,y) 



X\ X\ 



I 



+ I dxn{x) In 



n{x) 
TV- 1' 



(4.23) 



To find the n{x) which minimizes F, we take a functional 
derivative with respect to n and set it to 0, subject to the 
constraint that n{x) be normalized to A^ — 1. We obtain 
the following equation. 



dy V{x, y)n{y) -(- U{x) + V{xi , x) + In 



n{x) 

N - 1 



-A = 



(4.24) 
where A is the Lagrange multiplier enforcing the normal- 
ization constraint. We now turn to evaluating the eigen- 
value density, n{x). One can numerically solve this equa- 
tion, but it is possible to make reasonable approximations 
based on the qualitative picture that for x close to xi, 
the external potential determines the form of the density 
while for large x, the interaction dominates, crystalliz- 
ing the eigenvalues into a kind of lattice. Given the fact 
that we expect xi to be large, we may make the following 
approximation on the interaction. First, we may neglect 
the T^ interaction as it is rather small in comparison to 
the In sinh interaction. Secondly, we may make the ap- 
proximation that ln(sinh Xi — sinh Xj) « 2max(xi, Xj). 
On the external potential, we may make the following 
saddle point approximation: 



U{x) + V{x, xi) « C/,p(x) = U' + F'(x - 



(4.25) 



where, 



U' 

^ sp 

F' 



sp 



'— +ln2F 
4F 



F- 
1 

2f 



2F^ 

+ 1. 



(4.26) 



All of these approximations are quite good in the insulat- 
ing regime. Following this, our density equation becomes, 



7 / dyxn{y) +7 / dyyn{y) 



+t/,p(x)+ln^^+A = (4.27) 

This non-linear integral equation can now be reduced to 
a first order non-linear differential equation. To do so, 
we will first need the following boundary conditions im- 
posed by the integral equation itself. For x close to xi, 
the entropy (In) term will dominate over the interaction 
(integral) term and we'll have. 



i{x) ~ exp[-C/sp(a;)] 



Xi. 



(4.28) 



For large x, the interaction will dominate over the en- 
tropy, and we'll have: 



n{x) 



2r' 

7 



00 



(4.29) 



Now we take two derivatives of our integral equation for 
n{x). This results in the 2nd order non- linear differential 
equation: 

- n'ix) + -^P- - 2T'n{x) + -in{xf = 0. (4.30) 
n[x) 

Happily, this equation is autonomous in x, so we can 
change the independent variable to n, and the dependent 
variable to p = dn/dx. This will give a Bernoulli differ- 
ential equation in p. Using the boundary condition in 
Eq. (|4.29p . we can obtain the following differential equa- 
tion for n. 



dfi 

— =?i,/27n-4r'ln- , . 

dx V ' 2r'/e7 



(4.3f) 



The full solution to this equation is unknown, but in the 
insulating state, it is the small |a; — xij behavior that is 
important since only the lowest lying eigenvalues control 
the conductance. For x close to Xi, and therefore small n, 
the ln(n) term dominates the n term in the square root. 
Neglecting the linear term, we can perform the integral. 
Employing the boundary condition Eq. (|4.28p , we obtain. 



n{x) 



2V' 

67 



-T'(x^x'^^f 



X ^ Xi. 



(4.32) 



The numerical solution to this equation is shown in 
our previous paper— for two different values of F. It 
illustrates that the density is roughly constant with an 
exponential gap from the origin that increases with in- 
creasing disorder. In contrast, the density in the metallic 
regime is also constant, but starting at the origin. Thus 
our result suggests that the opening of a gap in the eigen- 
value spectrum could be considered as a signature of the 
metal-insulator Anderson transition. 



Now that we have determined the density of eigenval- 
ues for X close to 2:1, we can use it to evaluate Fp in this 
same region. We obtain, 



Fp{x,) 



8e7Vr' 



=erfc 



{xi - x')^ 



(4.33) 



Finally, this brings us to our result for p{xi). Using 
Eq. (|4.I8p . we obtain. 



p{xi) ex exp[-/(xi)] 



fix) 



U{x) 



(4.34) 



8e7Vr' 



zcrfc 



{x - x' )VF 



V. DISTRIBUTION OF CONDUCTANCE 

To obtain the probability distribution of the conduc- 
tance in the insulating state, we use the aformentioned 
fact that the dominant contribution to g comes from the 
first eigenvalue, xi. Therefore \Ti{g) ~ ln(f/cosh^(a;i)), 
which implies xi ~ (1/2) ln(4/g). And so, 



pilng) ex exp 



'2 g 



(5.1) 



This is our expression for the conductance distribution 
in the insulating state. In Fig. [SJ we compare the ana- 
lytical formula in Eq. (|5.1I) to numerical results and also 
to the distribution obtained from the GDMPK equation 
neglecting the interaction. 




Ing 



FIG. 9: (Color online) P{\ng) in the insulating regime for 
two values of disorder, F = 0.014 (solid red line) and F = 
0.045 (dotted green line), corresponding to i).ng) = —39.4 and 
(Ing) = —8.9, respectively. The numerical data points are 
from Ref. [2y ] for the same value of (In g) . The dashed blue 
curve is the GDMPK prediction of P{\ng) for (In 5) = -39.4 
neglecting the eigenvalue interaction 

Examining Eq. (j5.1|) we find the distribution has a 
Gaussian tail in Ing for very small g. We can calculate 
the expectation of In g as a function of F. The result, to 
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leading order, is 



(ln(5)> 



/4rin 



1 



4er3/2 



(5.2) 



One may also calculate the variance of lag in terms of 
(Ing). We find that asymptotically, 

(ln(i/g)) 



war (Ing) 



(5.3) 



ln(ln(l/g))- 

Current numerical simulations^ suggest rather that the 
variance of Ing goes as (— Ing)^'^. In our previous 
paper— we found a power law consistent with 2/3. It is a 
consequence of Eq. (15. 3p that different power laws can be 
fitted over different ranges of Ing. To illustrate we plot 
varlngvs. (— Ing)'^/"' inFig. (TTO]) . and in the insert, over 
a shorter range of Ing, we plot varlng vs. (— Ing)^'"^ 
which is actually a better fit over that shorter range. So 
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FIG. 11: The skewness is plotted vs. (—Ing). Note the 
asymptote at x ~ 1-1. 
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FIG. 10: Variance of Ing plotted as a function of (—Ing) ' . 
The points are calculated from Ea. (l5.1|l . The line is a best 
fit. Variance of Ing is plotted as a function of (— Ing)^'"' in 
the insert 

the discrepancy could be the result of the range of current 
numerical simulations (as the analytic solution seems to 
match with different power laws in different regimes), or 
the result of the series of approximations made in calcu- 
lating the partition function, which provides motivation 
for further work in this area. 

A major prediction of Eq. (|5.ip is the resulting novel 
disorder dependence of the skewness of the distribution. 
As shown in Ref. [2g, the skewness is positive in the 
deeply insulating regime and seems to tend to a constant 
value as disorder increases. From Eq. (|5.ip . one can cal- 
culate the skewness and confirm that it does approach 
a constant value, x(lii5) '^ l-l- The convergence to the 
asymptote is logarithmic, as illustrated in FigfTTl 

As disorder decreases the skewness eventually changes 
sign while still in the insulating regime. This change 
of sign occurs well before the Anderson metal-insulator 
transition point. Unfortunately due to the approxima- 
tions used to obtain the density of eigenvalues, we can 
not trust Eq. (j5.ip to describe the distribution at the 
transition point itself. 



VI. SUMMARY AND CONCLUSION 

In this work we start with the assumption, well sup- 
ported numerically, that the generalized DMPK equa- 
tion (GDMPK) with two scale-dependent parameters 
correctly describes the evolution of transmission eigenval- 
ues with increasing length in a 3D disordered conductor. 
Taking advantage of a mapping of the GDMPK equa- 
tion on to an imaginary time Schrodinger equation for 
a system of bosons interacting via a short range inter- 
action in ID, we then develop a systematic perturbation 
theory to solve for the full distribution of the eigenval- 
ues in terms of the A'^-particle Greens function describing 
the system (Eq. p. 141) ). In the strong disorder regime the 
strength of the interaction between bosons is small, which 
we use as the small parameter in our perturbation theory. 
In particular, we make sure that the TV-particle Greens 
function satisfies a highly non-trivial boundary condition. 
To leading order which includes at least all two-particle 
correlations, we obtain the distribution of the eigenval- 
ues in the strong disorder regime, where all eigenvalues 
are exponentially large. Writing the eigenvalue distri- 
bution in the form of a classical partition function, we 
then separate out the smallest eigenvalue x\ which dom- 
inates the distribution, and use a mean field approxima- 
tion that allows us to evaluate the density of the eigenval- 
ues (Eq. (|4.32p ) and the free energy close to x\. Finally, 
using Landauer formula connecting conductance with the 
transmission eigenvalues, we arrive at our expression for 
the full distriburtion of conductance, given by Eq. (15. ip . 
It turns out that the final distribution can be expressed 
in terms of a single disorder parameter. 

We use Eq. (|5.ip to compare directly with conduc- 
tance distributions evaluted numerically from three di- 
mensional tight binding Anderson model in the insulat- 
ing regime. Once the parameter of the theory is fixed to 
give the correct mean value of the distribution, the entire 
distribution is reproduced very well. In addition, we eval- 
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uate the variance and skewness as a function of disorder. 
While the numerical data for the variance is usually fitted 
with power laws, we find that the asymptotic behavior is 
perhaps more complex, given by Eq. (|5.3p . which can be 
fitted with different power laws in different regimes. We 
predict the skewness approaches a constant value in the 
deeply insulating regime, and changes sign as the disor- 
der is decreased from deeply insulating regime towards 
the Anderson transition point. Our expression becomes 
less applicable as one approaches the transition, but the 
change in sign occurs well before the transition. 

The theory can in principle be improved by keeping 
higher order terms in the perturbation expansion, as 
shown in Appendix A. However, even though the pertur- 
bation parameter remains small up to the Anderson tran- 
sition point, the partition function and hence the density 
of the eigenvalues can not be evaluated very accurately 
as the system approaches the transition. It is not clear at 
this point how to improve the calculation of the partition 
function beyond the mean field approximation used here. 
Nevertheless, starting at the strongly insulating regime, 
the present approach already predicts highly non-trivial 
changes in the conductance distribution as one decreases 
the disorder towards the Anderson transition. 



Appendix A: Corrections to Gn 

If we want to include three-particle correlations, then 
to leading order they would come from the following 
type of diagram (Fig. p2|) ) with two interaction lines. 
'Permutation' refers to the same graphs but with the or- 
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FIG. 12: tliree particle correlation diagram 

der of the first and second interactions switched. There 
would be the same set of diagrams for every distinct triple 
of particles within the set of N . The first dimensionless 
graph would correspond to an / of: 
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Let us call this term /3(1, 2, 3; 1', 2', 3'). Each of the other 
five graphs is simply a permutation of the arguments of 



G(kJoj) = 



ICD- £r 



FIG. 13: Unperturbed Green's function Go in Fourier space 



the first. For instance, switching the order of the in- 
teraction lines of /a would give us: /3(1', 2', 3'; 1, 2, 3). 
The second / would be /3(2, 1, 3; 2', 1', 3'), and its per- 
mutation would be /3(2',1',3';2,1,3). The third / 
would be 73(3, 2, 1; 3', 2', 1') and its permutation would 
be /3(3', 2', 1'; 3, 2, 1). Let i^3(l, 2, 3) denote the sum of 
these 6 /'s. Then using Eq. p.l4p we can write out the 
expression for Gn correct to second order in the interac- 
tion: 



Gi{i)y^G2{i,3) 



G--G^°^ns^n 



n 



G"(i)-liGO(i,j) 



:,F3{id,k) 



(A2) 



In order to calculate these diagrams it is convenient to 
go to a plane- wave basis. Unfortunately, the potentials in 
Eq. (|2.8I) do not possess well-defined Fourier transforms 
due to the l/x^ singularity at the origin. Even then, 
we have found that a good representation of the Fourier 
transform is given by the following relations: 



u{q) = Ai,7rg coth(7r(7/4) 
v{q) ~ \yT:qcoi}i{nq/2). 



(A3) 



The large q behavior of these potentials result in In- 
divergences of the first order (and probably higher order 
too) dimensionless diagrams, w. But the corresponding 
/'s result in well-defined, indeed necessary, terms by tak- 
ing the divergences to the exponent. This is one of the 
advantages of the exponential expansion, the divergences 
are resummed to give finite results order by order which 
we see accumulated in Gi and G2. On the other hand, 
if one is not interested in the small Xi or small Xi — Xj 
behavior (this ought to be taken care of via Gi and G2 
in any event), we may approximate these Fourier trans- 
forms by their small q limit 



u{q) 
v{q) f 



2A„ 



(A4) 



This is equivalent to approximating the potentials with 
a delta function which gives the same phase shift. Since 
the potentials are very short ranged, this should be ade- 
quate. Whichever set of expressions we use, the Feynman 
rules in frequency space would be given in Fig. (|13|) , and 
Fig. dUl). 

The first vertex corresponds to the interaction poten- 
tial v{q) — Xy/ sinh {xi — Xj), and the solid vertex corre- 
sponds to the interaction potential v{q) = A^/sinh (xi + 
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Xj). The latter can be incorporated into Fourier space 
provided we associate it with momentum q flowing from 
the center outwards towards both xi and X2, rather than 
with momentum q flowing from xi to X2 as with the dif- 
ference potential. Momentum would be conserved at all 
intersections as usual. 



FIG. 14: Interaction vertices in Fourier space 
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